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A new concept of using focus-diverse point spread functions (PSFs) for modal wavefront sensing (WFS) 
is explored. This is based on relatively straightforward image moment analysis of measured PSFs, which 
differentiates it from other focal-plane wavefront sensing techniques (FPWFS). The presented geometric 
analysis shows that the image moments are non-linear functions of wave aberration coefficients, but notes 
that focus-diversity (FD) essentially decouples the coefficients of interest from others, resulting in a set of 
linear equations whose solution corresponds to modal coefficient estimates. The presented proof-of-concept 
simulations suggest the potential of the concept in WFS with strongly aberrated high SNR objects in particular. 
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Optical wavefront is aberrated during its passage 
through physical systems due to perturbations like tur- 
bulence in Earth atmosphere or misalignment of tele- 
scopes. Useful understanding of the observed behaviours 
of the systems can be obtained by sensing the aberra- 
tion coefficients of the wavefront. Typically, a wavefront 
is measured at the pupil of an imaging system, as in 
Shack-Hartmann sensors (SHS) where wavefront slope is 
measured and then wave aberration is determined from 
the measurement pQ. It can be efficient and accurate, 
but requires a separete relay for pupil re-imaging with 
additional optical components. 

In some cases, one wishes to sense wavefront across 
the field of view of an instrument by using the built- 
in components without separate relay components. Such 
desires can be satisfied by applying FPWFS techniques 
to focus-diverse images. Curvature sensing (CS) is a FP- 
WFS method where the intensity difference between two 
extra- focal images are used to determine wavefront [2|. 
The CS signal is directly related to the shape of a mem- 
brane deformable mirror, enabling efficient control of the 
wavefront compensator pQ. In Phase retrieval, another 
FPWFS method, wavefront is optimized until its syn- 
thetic focus-diverse images closely match the measured 
ones [3]. It can be accurate and flexible, but may be less 
efficient than others due to a potentially large number- 
crunching during optimization. However, progress has 
been made to address this issue [U{6|. 

In this Letter, the idea of a new type of FPWFS tech- 
nique is explored. Like others, it utilizes focus-diverse 
PSFs, but differs, in that aberrations are sensed via im- 
age moment analysis of the PSFs. The presented geomet- 
ric analysis shows that the moments are non-linear func- 
tions of wave aberration coefficients, but notes that FD 
essentially decouples the coefficients of interest, leading 
to linear equations whose solution corresponds to modal 
coefficient estimates. It is argued that this approach has 
potential in WFS as shown in the presented simulations. 



Noll's notations are adopted throughout [7]. 

Geometrically, aberrated systems spread rays around 
and blur their PSFs, but the blur shape uniquely dif- 
fers depending on wave aberrations, e.g. comet-like PSF 
blur by coma (see pp. 90 in [8, for more examples). Such 
geometric relations can be explained by the connection 
between ray and wave aberrations. Suppose that a wave- 
front on the unit-disk pupil Q is given as, 

$ = ai Zi + a 2 Z 2 + a 3 Z 3 H h a M Z M , (1) 

where are the wave aberration coefficients and Zi are 
the Zernike polynomials in the pupil coordinates (x,y). 
The ray coordinates (X,Y) at the focus are given by 
differentiating by x and y, respectively. 

X = X + 2F— , F = y + 2F— , (2) 

where (Xq, Yq) are the ideal image coordinates and F is 
the focal ratio. The derivative terms are given as, 
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where a ,x and a ,y are slope coefficient vectors and Z is 
the vector of Z{. Noll noted the derivatives of Zi can be 
given in terms of Zi through the conversion matrices T x 
and T y [Jj. a!* and o!\ are related to a^, e.g. for M=10 
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One way to quantify the PSF shape is computing its 
kth. moment of order k = n + m > 1, that is given as, 
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where / is the pupil illumination and S is its integral 
over Q. Here, the moment is computed using the pupil 
plane quantities (i.e. J, X, and Y), whereas this is usually 
done using focal plane quantities. But, note that the focal 
plane PSF is essentially determined by the pupil plane 
quantities. Thus, Eq. [5] should lead to the moment as 
computed in the usual way. 

With 1 = 1 over fi, S = tt and the centroids become 
Mio = a'i and fioi = o!\. Apparently, ji nrn with k > 2 is 
non-linear in a'* or a'f , but only some slope coefficients 
(say bi) are coupled to {cL , 2)° t { a %Y wnere + ^ + 1 = k. 
Note that a' % and a'\ are the only coefficients affected 
by defocus (0,4), thus differentiating these with respect 
to <24 by k-1 times decouples bi from (a'^^i^Y- For 
instance, the derivatives of fi nm at k=2 (Eq. [6]), 
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from Eq.|U Another example is for k = 3, 



where the 2nd partial derivatives of /i nm are given as 
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where c = 12\/3 and d = 4\/3. Solving the above equa- 
tions produces the slope coefficients. In a similar way, 
the A:- 1th partial derivatives of the kth moments can be 
used to estimate the slope coefficients at the k- 1th or- 
der. The slope estimates can then be used in Eq. [3] to 
determine a$ of radial order up to k. 

What does it mean to have the derivative of moment 
in reality then? It is simply FD. Through- focusing is one 
way to achieve it, where the PSF is recorded at N differ- 
ent intra- or extra-focal planes and the image moments of 
the PSFs are measured. This obviously leads to /i nm that 
is a function of Aa 4 . Note N=k and Aa 4 =AL/(16V3F 2 ) 
in the unit of focus perturbation AL. Finally, fitting 
a polynomial of degree k-1 to the measured fi nm re- 
sults in the k-lth coefficient that corresponds to the k- 
lth derivative of /i nm with respect to 0,4. In summary, 
one can obtain the dependence of /i nm on a^, o!\ from 
Eql2l [3l H] and the fc-lth derivative of /i nm with respect 
to <24 leads to simple expressions allowing to estimate 



di. To experimentally measure the derivative, the PSF 
should be recorded at k different image planes or FDs. 

As a proof-of-concept, I demonstrate a simulation 
where non-zero with i = 2,3,4,---,10 are esti- 
mated by using three through-focus images sampled at 
-1mm, Omm, and 1mm (Aa4~±0.7A) from the ideal fo- 
cus of a f/10 system with 2.7m aperture (plate scale 
~ 7.6"/ram). A standard FFT method was used for syn- 
thesising polychromatic PSFs [9|. 




Fig. 1. Polychromatic through- focus PSFs on 512 2 grid with uni- 
form spectral weight at 11 wavelengths between 514nm and 614nm, 
Top: Initial, Bottom: After correction, Left: -1mm before focus, 
Middle: at focus, Right: +lmm after focus. 50/ira squares overlaid 
for size comparison. Images in square-root scale. 

The initial wavefront error is quite substantial 
amounting to 5. 7 A in rms at A=514nm, making the 
through-focusing effect in the PSF shape barely notice- 
able. The pixel size is 5/im. The estimated coefficients 
are compared to the input values in Table [TJ The esti- 
mated slope coefficients (a^, a!\ ) match the true values 
(a^, a!\ ) with sub-pixel accuracy. The maximum differ- 
ence between the estimated (di) and the true (a;) wave 
coefficients is 0.032 A in ag. The PSFs after applying di 
as corrections are in the bottom row of Figure [U 



Table 1. Slope (in pixel) and wave (in A) coefficients in Figure [T] 





a' X i 


«' x i 


a i 


IV 

a i 


di 


di 


l u 


-11.926 


-11.959 


-1.078 


-1.144 






2 1 


-27.885 


-27.986 


7.956 


7.968 


-1.198 


-1.185 


3 1 


7.956 


7.968 


-30.874 


-31.015 


1.612 


1.636 


4^ 


-12.351 


-12.505 


-13.033 


-13.316 


-4.241 


-4.241 


5 2 


1.823 


1.766 


-18.399 


-18.601 


1.624 


1.626 


6 2 


0.932 


0.917 


20.255 


20.598 


0.305 


0.309 


7^ 










-1.330 


-1.359 


8 3 










-1.261 


-1.276 


9 3 










1.593 


1.614 


10 3 










1.392 


1.408 



In the next demonstration, the same images were 
binned by a factor of 4 (Figure [2]). The same moment 
analysis was applied to the binned images and the resul- 
tant wave and slope aberration coefficients are compared 
to the true values in Table O The binning effectively re- 
duces the spatial resolution of the images. This leads 
to a fewer pixels used in the moment analysis and thus 
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reduces the moment measurement accuracy. To some ex- 
tent, binning is equivalent to low-pass filtering, washing 
out small-scale diffraction structures that may not be 
accounted for by the geometric wavefront aberrations, 
although diffraction effect appears to be less of an is- 
sue in the current analysis. Overall, the binning impact 
seems more prominent in the coefficients with odd radial 
orders (k=l, 3) than in the even radial order coefficients 
(fc=2). The estimated slope coefficients still approximate 
the true values with sub-pixel accuracy. 
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Fig. 3. Four object images at 1A FD (inset), RMS residual error 
against SNR for 21 objects. 




Fig. 2. Polychromatic through- focus PSFs with 4x4 binning, Top: 
Initial, Bottom: After correction, Left: -1mm before focus, Middle: 
at focus, Right: +lmm after focus. squares overlaid for com- 

parison. Images in square-root scale. 

The maximum difference between hi and a$ is 0.032 wave 
in clq. The through- focus PSFs after applying hi as cor- 
rections are shown in the bottom row of Figure El 



Table 2. Slope (in pixel) and wave (in A) coefficients in Figure 
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Finally, random aberrations (from a 2 to ai5, i.e. 
M=15) were applied to 21 point objects with different 
brightness. The aberration, without tip/tilt, amounts 
to 0.5A rms at 500nm. The images of four object are 
shown at the maximum FD of 1A in the inset (Fig- 
ure [3]). Five FD frames were obtained, where the pixel 
size is 18/ira with scale of 0.15" at 2.7m aperture and 
uniformly weighted 11 wavelengths between 500nm and 
600nm were used with photon shot noise. The signal- 
to- noise ratio (SNR), the square root of total photons 
collected as computed at the zero FD, ranges from 11 to 
221. The rms residual error, i.e. the root sum square of 
hi — di, is plotted against SNR (blue) with two power- 
law fits. While the green curve fits the blue data better, 



the red fit also asymptotes the data and follows the gen- 
eral behavior of the accuracy of other WFS techniques 
(~SNR _1 ) [l]. It appears that SNR higher than 50 per 
frame is needed to obtain accuracy better than 0.1 A. 
This translates to SNR-20 for a SHS with N=30 sub- 
apertures, based on the requirement of 2M<N pT| [T2]. 
leading to the theoretical SHS accuracy of 0.05A pQ. The 
analysis indicates the proposed technique to be suitable 
for high SNR cases, while spatial filtering, noise reduc- 
tion, or thresh-holding could help improve its robustness 
in low SNR cases. 

The demonstrations suggest the feasibility of the pro- 
posed geometric concept in WFS with high SNR and 
strong aberrations in particular. The fact that the mo- 
ment analysis is just an extension of the routine com- 
putation of image centroid or full-width-half-maximum 
can make all appropriate objects recorded in an image 
frame be potential WFS targets. This can certainly per- 
mit a straightfoward and rapid WFS potentially over a 
wide field of view, which can be particularly attractive in 
active diagnosis of alignment-driven field aberrations of 
imaging systems [13] or wide-field adaptive image com- 
pensation. The error propagation is to be separately dis- 
cussed with other considerations including aberration 
aliasing, diffraction, and FD implementation. 
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